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We investigate the dynamics of solitons of the cubic Nonlinear Schrodinger Equation (NLSE) 
with the following perturbations: non-parametric spatio-temporal driving of the form f(x,t) = 
a exp[i K (t) x], damping, and a linear term which serves to stabilize the driven soliton. Using the 
time evolution of norm, momentum and energy, or, alternatively, a Lagrangian approach, we develop 
a Collective-Coordinate- Theory which yields a set of ODEs for our four collective coordinates. These 
ODEs are solved analytically and numerically for the case of a constant, spatially periodic force 
f(x). The soliton position exhibits oscillations around a mean trajectory with constant velocity. 
This means that the soliton performs, on the average, a unidirectional motion although the spatial 
average of the force vanishes. The amplitude of the oscillations is much smaller than the period of 
f{ x )- I n order to find out for which regions the above solutions are stable, we calculate the time 
evolution of the soliton momentum P(t) and soliton velocity V(t): This is a parameter representation 
■ of a curve P(V) which is visited by the soliton while time evolves. Our conjecture is that the soliton 

becomes unstable, if this curve has a branch with negative slope. This conjecture is fully confirmed 
by our simulations for the perturbed NLSE. Moreover, this curve also yields a good estimate for the 
soliton lifetime: the soliton lives longer, the shorter the branch with negative slope is. 
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I. INTRODUCTION 



The Nonlinear Schrodinger Equation (NLSE) is one of the paradigms of soliton physics, because 
it represents a completely integrable system and has very many applications in practically all fields 
of physics, which are listed and discussed in several review articles [H, 0, 0] . For applications it is 
£2 ' important to study the perturbed NLSE 



iu t + u xx + 2\u\ u = R[u(x, t); x, t]. (1) 

Many different kinds of perturbations R have been considered and in particular the dynamics of a 
single soliton under these perturbations was investigated [3, 0]. 
In this paper we consider the following combination of perturbations 

•rH . 

X ■ R = f(x,t) -i/3u{x,t) - 6u(x,t) (2) 

with the real parameters (3 and 8 and the non-parametric spatio-temporal driving force 

f(x,t)=ae iK ^, (3) 

which yields several interesting effects, as we will see. The literature so far has mostly dealt with 
parametric driving [H, 0, 0, M, @, 0, HI]. Non-parametric (external) driving was studied without 
space dependence, e. f = eexp(iut) 0, Gil, EH, and with a periodic space dependence, e.g. 



/ = eexp[i(kx — cut)] |12|, ll3j]. Moreover, / = eexp[ig(x,t)_— iut], where g is a function of x — vt, 



was considered, but no localized solutions were discussed [13 . 

Our driving term in ([3]) was already used in the discrete form f n (t) = aexp[in<f)(t)]. Here the 
integer n denotes a lattice site in a nonlinear optical waveguide array which is modeled by a Discrete 



Nonlinear Schrddinger Equation (DNLSE) [14J. <f)(t) is the incident angle of a laser beam. In order 
to obtain a ratchet effect, a biharmonic <p(t) was used which breaks a temporal symmetry. 

In this paper we work with an arbitrary function K(t) in the driving term ([3]) and develop a 
Collective Coordinate (CC) Theory for the soliton dynamics which results in a set of nonlinear 
coupled ODEs for the CCs (Sections II and III). In order to obtain analytical and numerical solutions 
we then consider the case of temporally constant, spatially periodic driving f(x) = a exp (iKx) 
with constant K (Section IV). Although the spatial average of f(x) vanishes, there is transport: 
the soliton performs a unidirectional motion on the average, in contrast to the case of the driving 
R = V(x)u(x,t) with a periodic potential V(x) in which the soliton performs an oscillatory motion 
around a minimum of V(x) [7|. Solutions of the CC-Eqs. for the case of a harmonic or biharmonic 
time dependence of K(t) will be presented in a second paper. 

The second term, —if3u(x,t) with (3 > 0, in the perturbation (T2]) is a damping term which allows 
us to obtain a balance between the energy input from the driving and the dissipation. Other more 
complicated damping terms have been considered in [HI, EB]. 

The third term, —5u(x,t), in Eq. ([2]) will turn out to be decisive for the stability of the driven 
soliton. For 5 > the soliton radiates phonons (i.e., linear excitations) and eventually vanishes, 
or even breaks up into several solitons. For 5 < the situation is more complicated and will be 
discussed below. 

Section V presents a stable and an unstable stationary solution for the case without damping. 
For the region around the stable solution the CC-theory yields solutions in which all CCs exhibit 
oscillations with the same intrinsic frequency Q. However, tests by simulations, i.e. numerical 
solutions of the perturbed NLSE, reveal that the oscillatory solutions are stable only for certain 
regions of the initial conditions. These regions become broader when 8 is more negative (Section 
VI). 

We conjecture that the stability of any of these oscillatory solutions can be predicted by our CC- 
theory by calculating the curve P(V), where P(t) and V(t) are the momentum and the velocity of 
the soliton, respectively. This means that every point on this curve is visited during one period of 
the oscillatory solution. The conjecture is that the soliton will become unstable in a simulation, if 
the curve P(V) has a branch with negative slope, and this is confirmed by our simulations (Section 
VI). Interestingly, the curve P(V) not only predicts whether the soliton is unstable, but it also allows 
us to estimate the soliton lifetime: This time is longer, the shorter the branch with negative slope is 
(Section VI). 

The stability criteria for NLS-equations in the literature cannot be applied to our oscillatory 



solutions: The criterion of Vakhitov and Kolokolov [17|, [18( was established for stationary solutions, 
and the criterion of Barashenkov 0, Ji| for solitons travelling with constant velocity. In 0, 19] the 
slope of a curve P(V), different from ours, is considered and decides about the stability. However, 
here each soliton solution is represented by one point on this curve, i.e. the curve represents a family 
of solutions with different velocities. 

Finally we show in Section VII that the kinetic and canonical soliton momenta are identical, and 
we analytically calculate the soliton and phonon dispersion curves. 

II. TIME EVOLUTION OF NORM, MOMENTUM, AND ENERGY 

Multiplication of the perturbed NLSE (DO) by u* and subtraction of the complex conjugate NLSE 
multiplied by u yields 

| + | =i(fl .„_ fl „.) (4) 



with the density and current density 

q = \u\ 2 , j = i(u* x u — u*u x ). (5) 

By integration of Eq. (Tj0) over x, assuming decaying boundary conditions, we obtain the time 
evolution of the norm, 



N= / dx\u\ 2 , (6) 



which is 



N = i dx(R*u-Ru* 



+00 



= -2{3N + i I dx{f*u- fu*) , (7) 
J —00 

where the dot denotes the time-derivative and the terms with S have dropped out. 

Multiplication of Eq. (TjQ) by u*, addition of the complex conjugate equation multiplied by u x , and 
integration over x yields the time evolution of the momentum 



/+00 
dx{R*u x + Rul 
-00 



where P is defined as 



P=\ I <l.r{„„:.- ,/',/,). (9) 



Multiplication of Eq. (CQ) by w£, addition of the complex conjugate equation multiplied by u t , and 
integration yield the time evolution of the energy 



E = - I dx(R*u t + Ru* t ) , (10) 

where 

E= I dx[\u x \ 2 - \u\ A ] . (11) 



Interestingly, for (3 = and time independent force, i. e. R = f(x) — 5u, the r.h.s. of Eq. ( ITUl) 
can be written as a time derivative 

a r+00 

E = -— dx[f*u + fu* -5\u\ 2 }. (12) 



Thus the perturbed NLSE (0Q) without damping possesses a conserved quantity for arbitrary f(x 



"+00 

Titot 



E tot = / dx[\u x \ 2 -\u\ A -5\u\ 2 + f*u + fu*]. (13) 
J —00 

f(x) can be interpreted as a constant (external) force, in contrast to the case of the driving 
R = V(x)u, where V(x) can be understood as a potential (in which the solitons move). In this case 
the conserved quantity is [3] (see also (2p| ) 

r+00 

E£*= / dx[\u x \ 2 -\u\* + V(x)\u\ 2 ]. (14) 



In Sections II VI and IVl we will show that a soliton under a constant periodic force f(x) = aexp(iKx) 
performs, on the average, a unidirectional motion, although the average of / vanishes. 
When we include the damping {(3 > 0), the time evolution of the total energy is 

/+oo 
dx[2\u x \ 2 - 4|w| 4 - 25\u\ 2 + f*u + fu*} . (15) 
-oo 



III. COLLECTIVE COORDINATE THEORY 

The 1-soliton solution of the unperturbed NLSE reads [l| 

u{x, t) = 2ir] sech[2r/(x - £)] e -<0**+«) (16) 
with the real parameters rj > and £, the soliton position 

CW = Co-4^, (17) 

and the phase of the internal oscillation 

0(t) = 0o + 4(£ 2 -77 2 )t. (18) 

The soliton has amplitude 2rj, width l/(2rj), velocity V = — 4£ and phase velocity of the carrier wave 
V ph = -2£ + 2r) 2 /£. 
Including the term with 5 on the r.h.s. of Eq. (pQ), only the phase is changed: 

0(t)=0o + [4(e 2 -^ 2 )-5]t. (19) 

For this reason the term with 5 need not be treated as a perturbation and will be counted among 
the unperturbed parts of the NLSE in the following 

iu t + u xx + 2\u\ 2 + Su = f(x, t) — i(3u . (20) 

Therefore we now define the energy as 

/+oo 
cfe[H 2 -M 4 -%| 2 ], (21) 
-oo 

which is conserved for the unperturbed NLSE. Using Eq. f|T6|) we obtain the soliton energy 

16 

E so1 = 16^ 2 7] 3 - A5r] . (22) 

3 

The definitions ([6]) and (EJ) need not be changed, because all terms with S drop out on the r.h.s. 
of Eq. (J7J) and (jHJ). The norm and the momentum of the soliton are, respectively, 

N = Ar] ; P = -8r?£ . (23) 

Writing P = MV, the soliton mass is related to the norm by M = N/2. 

From the Inverse Scattering Theory (1ST) it is well known [lj that a perturbation theory results 
in a time dependence of rj and £ and in a change of the simple time dependence of ( and <fr i n Eqs. 
ffTTj) and ffTB"]) . Therefore we take Eq. f|T6l) as an ansatz with the four Collective Coordinates (CCs) 



r)(t),£(t), ((t) and <f>(t) (a similar ansatz has been used for optical solitons [2l|, |22j) and insert it into 
Eqs. ©, @ and $TD]), using f(x,t) = aexp[iK(t) ■ x] from Eq. ©. 



Eq. (JID yields 

f] = —2/3/7 - a^sechA • cos B, (24) 

with 

A(t) = ^[K(t) + 2m/v(t), (25) 

B(t) = <l>(t) + [K(t) + 2Z(t)]({t). (26) 

Two of the three terms which result on the r.h.s. of ([H]) cancel with the term — 8?)£ on the l.h.s.. 
The remaining terms are 

£ = aA sechA ■ cos B . (27) 

Finally, Eq. (FlUj) gives 8r]C,{...} — Ar]{...} = with two different curly brackets. This equation can 
be fulfilled by setting both curly brackets to zero, which yields 

C = -46 H -sechA • tanhA • sin 5, (28) 

877 2 

+ 2(6 = 4(£ 2 - rf) - 5 + ^— sechA • tanh A ■ sin B . (29) 

2?7 

We can alternatively use the Lagrangian method which will also be useful for section IVIII Our 
perturbed NLSE is equivalent to an Euler-Lagrange-Equation, generalized by a dissipative term, 

il DC ^ d dC dC dT 



dt ditt dx du* x du* du t 



with the Lagrangian density 



C = l -( Ut u* - u*u) - \u x \ 2 + M 4 + S\u\ 2 - fu* - f*u, (31) 
and the dissipation function 

T = i(3(uu* t - u*u t ) . (32) 

Inserting our CC-ansatz and integrating over the system we obtain the CC-Lagrangian 

16 

L = Ar](f) + 8r/C6 - 16r/^ 2 H r/ 3 + 45rj - 2vra sechA ■ sin B (33) 

3 

and the CC-dissipation function 

F = -P(8rj(j> + I677CO • ( 34 ) 
The CC-equations are then obtained by the four generalized Lagrange equations 

d dL _dL _ OF 

where ip stands for the four CCs rj, 6, C an d 4>- The resulting four ODEs are identical with the Eqs. 

(El, (EZD, (EHD and (EHD. 

Finally we evaluate Eq. ( fl~3j) . yielding 

16 

E tot = lQr]£ 2 T] 3 - A5r] + 2ira sechA ■ sin B . (36) 

3 

Here the first three terms are the soliton energy (1221) . while the last term stems from the perturba- 
tions. We note that E tot is conserved only in the case of no damping (/3 = 0) and time independent 
force f(x) = aexp(iKx) with constant K, see below Eq. (fT2l . 



IV. CONSTANT, SPATIALLY PERIODIC FORCE 



We consider this case because it exhibits some surprising, counter-intuitive features: E. g., the 
soliton position performs oscillations on a length scale that is very different from the spatial period 
L of the force. 

We take a constant K in Eq. ([3]), i. e. f(x) = aexp(iKx), and consider only small values of \K\ 
such that the period L = 2n/\K\ is much larger than the soliton width. We first consider the case 
with damping (/3 > 0) for which we can expect "steady-state solutions" (see below Eq. (HOI) ) for 
times much larger than a transient time r on the order of 1/(3. 

The transformation u(x,t) = ty(X,t) exp(iKx) into a moving frame X = x — Vft with Vf — 2 K 
leads to the autonomous equation 

i% + * XI + 2|*| 2 * = a + (K 2 - 6)V - (37) 

with the non-parametric constant driving term a. Apart from the factor K 2 — 5 := c 2 , which can 
be eliminated by scaling time by c 2 , space by c > 0, and * by 1/c, Eq. (]3~T1) is the same as an 
autonomous equation, which was obtained from the NLSE (JTJ) with R = eexp(iut) — i(3u by the 
substitution u = $?(x, t) exp(iut), setting u — 1 0, llll ]. However, these investigations differ from 
ours in several respects: In Ref. [Io| static soliton solutions of Eq. (1371) were obtained numerically, 
then an existence and stability chart was constructed on the (a, 0) plane. In Ref. [ll[ a singular 
perturbation expansion was performed at the soliton's existence threshold. In contrast, we study 
moving solitons by solving our CC-equations, and test the results by simulations, i.e. by numerically 
solving the NLSE (CQ). Besides the steady-state solutions for (3 > (this section), for the case 
(3 = we study stationary solutions and oscillatory solutions, which have a more complicated time 
dependence (Section V). 

The CC-equations from the previous section yield steady-state solutions, in which the driving is 
compensated by the damping, using the ansatz C = Co + Vft with Vf = 2K and constant rj, £, and 0: 

uf(x,t) = 2iri f sech{2rif[x - C(t)]}e i{Kx '^t \ (38) 

with 

rjf = VK 2 - 5/2, (39) 

$ = ±- + arcsin (40) 
1 2 na 

We denote this as steady-state solutions, because (f> is constant, in contrast to stationary solutions, 
where has a linear time dependence (Section V). These solitons have an internal structure due to 
the factor exp(iKx), but no internal oscillations since 0j is constant. In the moving frame these 
solitons correspond to the above mentioned static solutions of Eq. (}3"7]) in Ref. (Io| . 

We have numerically solved the CC-equations for many sets of initial conditions (IC) 770, Co ; £0 an d 
4>q. There is a basin of attraction around the solution (|38[) with <p^. The solitons always evolve 
to this solution, except when the values of r]o, £0 and 0o are too far from those of this steady-state 
solution and when the damping (3 is too large. An example for this is the parameter set a = 0.05, 
K = 0.01, 5 = -3, (3 > (3 C w 0.003 with the IC r] = 1, £ = Co = 4>o = 0. Here the soliton 
vanishes, i.e. its amplitude and energy go to zero while its width goes to infinity. In order to achieve 
a convergence to the stable steady-state solution one can either reduce (3 below the critical value (3 C 
(which depends on the other parameters and the IC), or go closer to rjf = 0.8660 by reducing rj by 
0.1, for instance, or go closer to 0t by choosing O = n/2. 



In this context it is interesting to consider the total energ y (f5oT) and its time derivative in which 
the CC-equations can be inserted. Using MATHEMATICA [23j we obtain 

E tot (t) = -/3{8r/[4(£ 2 - T] 2 ) -5]+ An A sechA ■ tanh A ■ sin B} . (41) 

This is indeed zero for the steady-state solution because £/ = —K/2 and thus A = 0. On the path 
to that solution, Eq. (jUJ) alternately exhibits both signs. I. e., the total energy performs oscillations 
which become smaller and smaller while approaching the final value Ej ot . 

According to Eq. ( |39i) . 5 < K 2 . As we choose \K\ -C 1 (see above), 5 is either positive but very 
small, or 5 is negative. This results from the CC-theory; in section IVlIl we will show that the driven 
soliton can be stable only for 5 < 0, by taking into account the phonon modes. 

The quality of the CC-theory must be tested by simulations. The soliton shape agrees very well 
(Fig. [T]), but in the simulations the soliton resides on a small constant background because the 
perturbation f(x) = aexp(iKx) does not vanish far away from the soliton. This background is 

u bg = -—e iKx , (42) 

with u>k = K 2 — 5 — i(3; see also the last term in Eq. 0771) . Plotting the real and imaginary parts of u, 
the spatial period L = 2it/\K\ is observed. The norm density \u\ 2 forms a shelf on which the soliton 
moves; the shelf height is quantitatively confirmed. The dynamics of the soliton is practically not 
affected by the background: the time evolution of the soliton position is identical in the CC-theory 
and the simulations (Fig. [2} right panel); only the soliton amplitude differs a little (left panel). 




200 



FIG. 1: Left panel: Soliton moving to the left for t* = 250,500. Simulations of NLSE (solid lines) and numerical solutions of 
CC-equations (dotted lines). Right panel: Real (solid line) and imaginary (dashed line) parts of u(x,t) for t — 500. Parameters: 
K = -0.1, a = 0.05, 5 = -1, /3 = 0.05, with IC Co = 0, Co = 0, 4>o = 1.69 and rjo = 0.5. 

All CCs exhibit decreasing oscillations with an intrinsic frequency Q which will be discussed in the 
next section. However, oscillations are not visible in the soliton position ((t) in Fig. [5] because the 
linear term dominates the time evolution. By reducing the damping (3 and by choosing a smaller 
time scale one can see the oscillations also in ((t), see Fig. El 

In any case, the soliton performs, on the average, a unidirectional motion although the spatial 
average of the periodic force f(x) = a exp(iKx) is zero. Thus, this is a ratchet-like system in which 
the translational symmetry is broken by the inhomogeneity f(x) in the NLS-equation. 

The period L = 2tt/\K\ of the inhomogeneity f(x) is not reflected in the soliton dynamics, because 
our Eqs. ([I])-® could be reduced to the autonomous Eq. (I37j) . For early times (t < r = 0(l//5)) 
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FIG. 2: (Color online). The amplitude and position of the soliton obtained from a simulation of the NLSE (red dashed lines) 
and from the numerical solution of the CC equations (solid lines). Parameters: K — —0.01, a — 0.05, 8 — —3, (3 = 0.01, with 
IC Co = 0, Co = 0, <fro = tt/2 and 770 = 1. 
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FIG. 3: (Color online). The amplitude and position of the soliton obtained from a simulation of the NLSE (red dashed lines) 
and from the numerical solution of the CC equations (solid lines). Parameters: K = 0.1, a = 0.05, S = —3, (3 = 0.001, with IC 
Co = 0, Co = 0, 0o = and 770 = VK 2 - 5/2. 

the soliton performs the above mentioned oscillations. However, the amplitude of these oscillations 
is much smaller than L (e.g., L = 62.8 for K = 0.1). When the strength a of the inhomogeneity is 
strongly increased (a = 1, for instance), ((t) exhibits a staircase structure. The height of the steps 
is larger than the amplitude of the above oscillations, but still much smaller than L. However, such 
values of a represent a strong perturbation for which the CC-theory can no longer be valid. Indeed, 
in the simulations the soliton soon becomes unstable. 



OSCILLATORY SOLUTIONS 



In order to study in more detail the intrinsic oscillations in the CCs (see previous section), we 
consider the case without damping {(3 = 0). Here oscillatory solutions are possible because the total 
energy (1361) is conserved, see Eq. (HIT) . This means that in Eq. (1361) oscillations of the soliton energy, 
Eq. (122]) . are compensated by the oscillations of the term 27rasechv4(t) ■ sin B(t) stemming from the 
perturbations. This is confirmed by inserting numerical solutions of the CC-equations into Eq. 



In order to obtain analytical solutions our approach is to look first for stationary solutions and 
then to consider small oscillations around them. For the stationary solutions we make the ansatz 
C = C + V s t, r] = r] s ,£ = £ s and = <p s - a s t. The CC-equations (|23) ; (|27|), <|25j>, (|25|) yield 



(ITT 

= — — sechA, ■ cos 5, (43) 
= aA-sechA, ■ cos 5, (44) 

CUT 2 

V s = -4£ s H - sechA s • tanh A s ■ sin B, (45) 

Srjf 

-a s = 4(£„ 2 - 7] 2 S ) -5 + sechA • tanh A s ■ sin B, (46) 



with 



A s = 1(K + 2Q/ Vs , (47) 
B{t) = <f) s + (K + 2Q( S + [(K + 2QV S - a s }t . (48) 
Eq. ( )43l) must be fulfilled for arbitrary t, thus cosS = 0, which leads to 

a s = (K + 2£ S )V S} (49) 
<f>. + (K + 2Z.)( a = ±|, (50) 
sin B = ±1, (51) 

respectively We distinguish two different cases: In case I, a s — and we obtain the steady- 
state solutions of Section IV. In the co-moving frame these solutions correspond to two exact static 
solutions of ( 1371) for zero damping I25H. Here one solution is stable below a critical driving; 



strength, whereas the other one is always unstable. 

In case II, a s 7^ 0, which means that we obtain stationary solutions. Here we can restrict ourselves 
to the case ( s = £ s = 0, because other values yield qualitatively similar results. Using Eqs. (|45|) . 
(l4~6"j) . (j4*7j) and (14^1) one transcendental equation for rj s remains: 

7i 2 aK 

Aril = -8 ± ^-sechA s • tanhA s . (52) 

4r) 2 s 

For the parameter set a = 0.05, K = 0.1 and 5 = —3 we get r] s = 0.866239 for <p s = tt/2 and 0.865811 
for cf) s = — it/ 2. The numerical solution of the CC-equations for these two cases reveals that in the 
former case we have a stable stationary solution, whereas in the latter case the solution is unstable. 

For the region around the stable stationary solution we expect that all CCs exhibit oscillations and 
we assume that these oscillations are harmonic if the amplitudes are sufficiently small. We choose 
0o — tt/2 and make the ansatz 

((t) = Vt~a c sin(fit), (53) 
V(t) = Vo + ar, [1 - cos(ftt)], (54) 
f(t) = -o e [l-cos(fi*)], (55) 
4>{t) = 4>o — at + <i0 sin(fit). (56) 

This ansatz takes into account that the soliton first starts its oscillatory motion (i.e., ( — Vt is linear 
in t for small t), and then the soliton shape changes (i.e., r] is quadratic for small t). 



Since we have considered the stationary solution that belongs to £ s = ( s = 0, we can neglect £ 
compared to K and obtain for B in Eq. (121)]) : 

B = 0o + (tty - A"a c ) sin(fit) + (liu - a)*. (57) 

Ka^ is one order smaller than a^, because \K\ <C 1, and we get 

B = <f) + a lf> sin(Ot) + (.KV - a)t. (58) 

Here we can distinguish two limiting cases: In case 1, the linear terms cancel and we have a pure 
oscillatory behavior. In case 2, the linear term dominates the oscillatory term which can then be 
neglected. Indeed, taking the second derivative of with respect to time and using (124j) . ( 1271) . 
( 1281) and ( 1291) one obtains = ci sin(S) + C2sin(2£>) + C3Cos(£>) + C4Cos(2S) + C5, where q with 
z = 1, 2, 5 are functions of the CCs. Considering the above approximations, and in addition when 
the terms of order of aKV(4r)Q + S)t are negligible, one realizes that ci, C2, C4 and C5 are small 
compared with C3 = Q 2 , 

tt 2 « 4vrar/ , (59) 

and hence, the Eq. for 0(i) reads 

0-ft 2 cos(0) = 0. (60) 

Integrating this equation once we obtain 

= ±ttV2 y/sm((f>) + C, (61) 

where C = 4>l/(2Vt 2 ) - sin(0 o ) « ((4^ + <5) 2 )/(87tgw7 ) - 1- Then, if C > 1, the constant term 
dominates and so <p(t) goes linearly, i.e. 0(t) = 7r/2 ± (4 t]q + 5) t. Notice that C ^> 1 implies 
(4r/ 2 + 5) 2 ^> 167rar7o, which can be solved for r/ , yielding for a = 0.05, 6 = — 3 and O = 7r/2, 770 <^ 
0.655 or 770 ^> 1.077. This is confirmed by numerical solutions of the CC-equations which exhibit 
an oscillatory behavior of <p(t) for 770 G [0.65; 1.07], and a linear behavior (plus small oscillations) 
outside of this interval. 

Finally we remark that in case II both V s in the stationary solutions and V in the oscillatory 
solutions differ from the value IK in case I. A transformation to a frame moving with V ^ 2K 
would not simplify the NLSE because terms with *f?x would appear in Eq. (13T|) . 



VI. STABILITY OF OSCILLATORY SOLUTIONS 



Our simulations reveal that the driven undamped soliton is stable only for a part of the set of 
solutions obtained by the CC-theory. Naturally we would like to predict, by using the CC-theory, 
which solutions are unstable and to understand what causes the instability. (The latter point will 
be discussed in the next section). 

For the NLSE with a general local nonlinearity, but without driving and damping, the Vakhitov- 
Kolokolov stability criterion [13, 18[ states that solitons are stable if dN/dA > 0. Here N is the norm 



and A the so-called spectral parameter in stationary solutions of the form u(x,t) = $?(x) exp(iAi). 
However, our oscillatory solutions, Eqs. (|53l -( l56l) inserted into Eq. (fT6l) . have a more complicated 
time dependence than stationary solutions; therefore the criterion cannot be applied here. 

The same holds for the stability criteria of Barashenkov 0, which were established for solitons 
travelling with constant velocity. But here we can get some motivation for how to proceed in our 
case: Barashenkov showed that dark solitons of the NLSE with generalized nonlinearity are stable if 



dP/dV < (here the tildes are used to distinguish from our P and V) [19]]. Thisproof was carried 
over to (bright and dark) solitons of the undamped parametrically driven NLSE [6|. Here the point 
dP/dV = separates a stable from an unstable branch of the curve P(V), but it depends on the 
type of the solution on which side the stable branch is. For the following it is important to note 
that this curve P(V) represents a family of solutions with different velocities, i.e. each solution is 
represented by one point on the curve. 

We make the conjecture that our oscillatory solutions are dynamically unstable, if our curve P(V) 
has a branch with negative slope, i.e. 

dP , . 

JV < (62 > 

for a finite interval of V. This curve is obtained from its parameter representation P(t), V(t), where 
P(t) = —8r]^ is the soliton momentum Eq. ( 1231) and the soliton velocity V(t) = ( is obtained by the 
r.h.s. of Eq. (12811 . Each oscillatory solution is represented by its own curve P(V). This curve has 
a finite length, because P(t) and V(t) are periodic and remain finite, (C(i) and 4>(t), which contain 
terms linear in t, do not appear in P(t), and in V(t) they appear only via sinS, see Eq. (J28J). 
Plotting the "stability curve" P(V) we can inmediately see whether there is a branch with negative 
slope. 

For the parameter set a = 0.05, K = 0.1, (3 = and 5 = — 1 and IC Co = £o = and (fio = ir/2, 
we find a small "stability interval" 0.48 < rjo < 0.52, i.e. an interval of initial conditions for which 
the solutions are stable. As expected, this interval is situated around the value rj s = 0.501874 from 
the stable stationary solution (1521) . There is another stable regime for rjo > 0.76. When we go far 
away from the IC for the stationary solution by choosing 0o = 0, instead of <fio = n/2, the stability 
interval around rj s vanishes. The upper stability regime exists now for r) > 0.69. 

When \6\ is increased, e.g. by choosing 5 = —3, the stability interval around rj s = 0.866239 is 
much larger than in the case 5 = — 1; for <p = tt/2 it is 0.7 < t] < 1.03. The upper stability 
region is above rj = 1.08. Considering again O = 0, the stability interval around 7] s only shrinks to 
0.76 < r] < 0.97, but does not vanish because it was much larger than in the case 5 = — 1. The upper 
stable region is above t]q = 1.02. Thus the conclusion is that an increase of \5\ widens the regions 
with stable soliton solutions. All the stability regions given above are confirmed by simulations for 
the perturbed NLSE, with an error of less than 1 %. 

At some of the boundaries of the above stability intervals there is a drastic change in the shape of 
the solutions of the CC-equations and the stability curve. E.g., if we choose r] = 0.75 (Fig. HJ), which 
is just below the stability regime 7] > 0.76 (see above), we obtain very anharmonic oscillations in all 
CCs and the stability curve has a long branch with negative slope (Fig. The simulations indeed 
show that the soliton becomes unstable very quickly and vanishes (Fig. HJ). However, a slight change 
of i] to the value 0.76 produces a stability curve that has only one branch with a positive slope (Fig. 
[5]). Here the soliton is indeed stable in the simulations (Fig. [6]). 

At the boundaries of other stability intervals the change of the shape of the solutions is less 
dramatic. E.g., for the case 5 = — 1 and rjo = 0.46 the oscillations in the CCs are small and 
harmonic. The stability curve has only a very short branch with negative slope which is visited in 
the time evolution only for very short time intervals (Fig. [7j) . Here the soliton is indeed stable for 
a relatively long time (Fig. [H]). This soliton lifetime is increasingly reduced when 7] is reduced by 
which the negative slope branch in P(V) becomes longer. Thus this curve predicts not only whether 
the soliton is stable, but also gives an estimate for its lifetime when it is unstable. 




FIG. 5: "Stability curve" P versus V. Left panel: rjo = 0.75. Right panel 770 = 0.76. Other parameters as in Fig. [3] 
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FIG. 6: The evolution of rj, £, £ and cf> obtained from a simulation of the NLSE (dashed lines) and from the numerical solutions 
of the CC equations (solid lines). Parameters: K = 0.1, a = 0.05, 5 — — 1, (3 — 0, with IC £o ~ 0, Co ~ 0, <f)o = t/2 and 
J70 = 0.76. The simulations were carried forward to a final time tf = 1000. 



VII. SOLITON AND PHONON DISPERSION CURVES 



For the soliton dispersion curve E sol (P) we need the canonical momentum P of the soliton. P 
must fulfill the Hamilton equation 

C = Qp, (63) 
where ( is the soliton position. We start from the CC-Lagrangian Eq. ( 1331) and obtain 

It = ^ = N > ( 64 ) 

d(j) 

as the angular momentum of the internal oscillation of the soliton, and 

dL 

= 8 V ( = 2N( := X, (65) 

as the variable which is canonically conjugate to £. Using a Legendre transformation we obtain the 
Hamilton function 



H = X£ + N<p-L, 
H(X, £; N, 0) = AN£ 2 - ^A^ 3 - 5N + 2vra sechA • sin B, 



(66) 
(67) 
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FIG. 7: "Stability curve" close to the boundary of an unstable regime. Parameters: K = 0.1, a — 0.05, S — —1, /3 = 0, with 
IC £ = 0, Co = 0, 0o = tv/2 and r? = 0.46. 
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FIG. 8: The evolution of rj, f , £ and obtained from a simulation of the NLSE (dashed lines) and from the numerical solutions 
of the CC equations (solid lines). Parameters as in Fig. [7] 

with 



(68) 



We now make an ansatz for a canonical transformation to the following new set of variables: 



P = -2iV£, C = ^ X ( 69 ) 
N = N, = + s(iV,£,X). (70) 



This means that P is identical to the kinetic momentum in Eq. ( 1231) . N is chosen to be equal to N 
and g still has to be determined such that the four fundamental Poisson brackets are fulfilled. This 
yields g = £X/N '. The new Hamiltonian reads 

H(P, C; N, 4>) = ^P 2 - -^N 3 -5N + 2na sechA ■ sin B, (71) 



A = T-r ( K - Tf ) > B = 4> + K(. (72) 



with 



N \ N 

The first two terms in Eq. ( J7T1) agree with literature results on the unperturbed NLSE |3|]. We have 
checked that the four Hamiltonian Eqs. which result from Eq. ( 1TTT) are indeed equivalent to the 
four CC-equations in section II III for the case without damping. The above results hold for the force 
f(x,t) in Eq. with arbitrary K(t). In the following we return to the case of constant K in order 
to calculate analytically the stability and dispersion curves wherever it is possible. 

Let us consider the stability intervals around the stationary solutions as given in Section IVT1 
Except for the close vicinity to the boundaries of the stability intervals, the oscillations in the CCs 
are nearly harmonic and can be well approximated by the Eqs. ( 1531) to ( 1561) . Inserting into P = —8 r] £, 
neglecting the 2nd harmonic, and using V — £, one can easily see that P(V) is a straight line with 
slope 

because for r] < r] s both ag and are positive and for r] > r] s both are negative. The latter also 
holds for the upper stability intervals r]o > 0.76 for 5 = — 1 and 770 > 1.08 for 5 = —3. Near or at the 
boundaries of the stability intervals, the oscillations in the CCs are very anharmonic and therefore 
the calculation leading to Eq. ( 1731) is not possible. P(V) is a curved line which can be calculated 
using the numerical solution of the CC-equations. 

We now turn to the soliton dispersion curve E sol (P). Both P = — 877^ and E so1 in Eq. (1221) 
consist of powers of r\ and £. For simplicity we write Eqs. (j31D and (135]) as r\ = fj — a v cos(Ot) 
and £ = £ + cos(Qt), where £ is negligible. We concentrate on the leading terms in cos(Ot) and 
distinguish two cases: In case I, both, E sol (t) and P(t) have a leading term with cos(f2t); in this 
case the dispersion curve E sol (P) is linear in the first approximation. In case II, the l s *-order terms 
in E sol {t) cancel if fj = | V—S, but cannot cancel in P(t). In the next order 

E so \t) = E max - AE cos 2 (fit), (74) 

with AE = 16fj(a 2 /3 — a|) > because \a^\ \a v \. Inserting P = P + Sfja^ cos(Ot) yields a 
parabolic dispersion curve 

E so \t) = E max - AE (^r^) ■ (75) 



This turns out to be a surprisingly good approximation when comparing with the dispersion curve 
obtained by using numerical solutions for r](t) and even when the cancellation of the linear 
terms in E sol (t) is not exact. The condition fj = \J—8/2 is approximately fulfilled for the regions 
with strong oscillatory terms in see end of Section IVl 

When solitons become unstable they radiate phonons (i.e., linear excitations). Therefore we con- 
sider the perturbed linearized NLSE without damping 

iu t + u xx + 8u — a e lKx , (76) 

which is solved by 

u(x, t) = c e ^™*) + b e^*-"**) - — e iKx (77) 

with u k = k 2 — 5, u K = K 2 — 5. The first term in Eq. (ITT)) represents the phonons of the unperturbed 
equation with the free amplitude c and the dispersion curve u^. 

The second term in Eq. (I77j) with the free amplitude b represents a single phonon mode whose 
wave number and frequency are given by the parameter K in the force f(x) = ae tKx . Such phonons 
are typically radiated at the beginning of a simulation as the initial soliton profile adapts to the 
system. These phonons can be observed best when they interact with the soliton after having been 
reflected by a boundary of the system. 

Finally, the last term in Eq. ([77]) represents a static background with a fixed amplitude a/u K . 
This was already discussed below Eq. 

VIII. CONCLUSIONS 

We have considered the dynamics of NLS-solitons in one spatial dimension under the influence of 
non-parametric spatio-temporal forces of the form f(x,t) = a exp[i K(t) x], plus a damping term 
and a linear term 8 u(x, t) which stabilizes the driven soliton. We have developed a CC-theory which 
yields a set of ODEs for the four CCs (position (, velocity £, amplitude rj, and phase 0). 

These coupled ODEs have been solved analytically and numerically for the case of a constant, 
spatially periodic force f(x) = a exp[i K x]. The soliton position exhibits oscillations around a mean 
trajectory £ = Vt; this means that the soliton performs, on the average, a unidirectional motion 
although the spatial average of the force vanishes. The amplitude of the oscillations is much smaller 
than the spatial period L = 2n/\K\ of the inhomogeneity f(x). The other three CCs also exhibit 
oscillations with the same frequency as ({t). 

In the case of damping, the above oscillations are damped and the solution approaches a steady- 
state solution with constant velocity, if the IC are close enough to those of the steady-state solution 
and if the damping is not too large. Otherwise the soliton vanishes, i.e. its amplitude and energy 
go to zero while its width goes to infinity. 

In the case without damping all the above oscillations persist. These periodic solutions exist 
because the total energy of the perturbed system is a conserved quantity, even for arbitrary inhomo- 
geneity f(x), and independent of the CC-ansatz. 

However, a comparison with simulation results for the perturbed NLSE reveals that only part of 
the above oscillatory solutions are stable. Our CC-theory predicts the unstable regions in the IC 
and the parameter 8 with high accuracy, by using our conjecture that the soliton becomes unstable 
if the slope of the curve P(V) becomes negative somewhere: here P(t) and V(t) are the soliton 
momentum and velocity, respectively. It turns out that the stability intervals become broader when 
the parameter 8 is chosen more negative. Moreover, we have found that the curve P{V) also yields a 



good estimate for the soliton lifetime: The soliton lives longer, the shorter the negative-slope branch 
is, as compared to the length of the positive-slope branch. 

Other cases of the force f(x,t) = exp[i K(t) x] will be considered in a second paper: specifically 
single and biharmonic K(t), with and without damping. 
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